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Abstract 

Two-photon and two-gluon annihilation rates of r] c , rj r c , r]b and n' b are estimated in the relativistic 
Salpeter method. By solving the full Salpeter equation with a well defined relativistic wave function, 
we estimate M Vc = 2.979 ± 0.432 GeV, = 3.566 ± 0.437 GeV, M Vb = 9.364 ± 1.120 GeV and 
= 9.941 ±1.112 GeV. We calculated the transition amplitude using the Mandelstam formalism 
and estimate the decay widths: T(n c -► 2 7 ) = 7.14 ± 0.95 KeV, T(n' c -► 2 7 ) = 4.44 ± 0.48 KeV, 
T(n b -> 27) = 0.384 ± 0.047 KeV and T(n' b -» 2 7 ) = 0.191 ± 0.025 KeV. We also give estimates of 
total widths by the two-gluon decay rates: T tot (r] c ) = 19.6 ± 2.6 MeV, T tot (n' c ) = 12.1 ± 1.3 MeV, 
T t at(rib) = 6.98 ± 0.85 MeV and T tot (rj b ) = 3.47 ± 0.45 MeV. 
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I. INTRODUCTION 



It is well known that two-photon or two-gluon annihilation rate of heavy h quarkonium 
cc or bb is related to the wave function, so this process will be helpful to understand the 
formalism of inter-quark interactions, and can be a sensitive test of the potential model. 
With the replacement of the photons by gluons, the finial state becomes two gluon state, 
which will be helpful to give information on the total width of the corresponding quarkonium. 

Experimentally there are quite many results for the decay width T to t{Vc) with a wide range 
of values and uncertainties by different collaborations; for example, in recent experiment of 
Barbar Q] they give 34.3 (2.3) (0.9) MeV, much larger than the cited value 16.0iH MeV by 
Particle Data Group 0]. However, r]' c has been just declared observed by Belle jy] and by 
Barbar [1]; r/f, and r]' b have not been observed yet, even though there were some experiments 
to search for e.g., the ALEPH j^J collaboration. In short, unlike the corresponding vector 
1 quarkonium which can be produced directly by e + e~ annihilation, experiments on _+ 
quarkonium have just begun, even for i] c . This due to the small cross section; presently 
there are 57.7 x 10 6 J/\l/ events collected with the BES-II detector |5[, but there are only 
2547 ± 90 7] c events collected by the Barbar detector Q| . 

For the theoretical estimates of the annihilation rate for rjp and r)l\ we have various 
methods readily available in hand. First was the non-relativistic calculation, then the rela- 
tivistic corrections were found to be important especially for cc states. In recent years, many 
authors try to focus on the relativistic corrections and there are already some versions of 
relativistic calculation, and they give improved results over the non-relativistic methods. In 
this letter, we give yet another relativistic calculation by the instantaneous Bethe-Salpeter 
method |6], which is a full relativistic method [7] with a well defined relativistic form of the 
wave function. 

There are two sources of relativistic corrections; one is the correction in relativistic kine- 
matics which appears in the decay amplitudes through a well defined form of relativistic wave 
function (i.e. not merely through the wave function at origin); the other relativistic correc- 
tion comes via the relativistic inter-quark dynamics, which requires not only a well defined 
relativistic wave function but also a good relativistic formalism to describe the interactions 
among quarks. 

The Bethe-Salpeter equation is a well-known tool to describe a relativistic bound state. 
And the Salpeter equation is the special case of Bethe-Salpeter equation when the interaction 
is instantaneous. It has been shown that the instantaneous approach is a good approximation 
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in heavy mesons, especially for the equal-mass quarkonium, since the non-instantaneous 
correction was found to be very small in equal mass system j8j. The full Salpeter equation 
includes two parts of the wave function, the positive and negative energy part. In the case 
of heavy mesons the negative energy part usually gives a smaller contribution than the 
positive energy part, and therefore, to simplify the calculation for the heavy mesons authors 
like to make a further approximation to the Salpeter equation by ignoring the negative part 
contribution. However, since we are considering full relativistic calculation, and the negative 
energy contribution was found to be not very small for some cases jsj, in this letter we will 
solve the full Salpeter equation including the negative contribution, and use the full Salpeter 
wave function to estimate the annihilation decay width of quarkonium. 

We note that the form of the wave function is also important in the calculation, since 
the corrections of the relativistic kinetics come mainly through it. We begin from the 
quantum field theory, analyze the parity and charge conjugation of bound state, and give a 
formula for the wave function that is in a relativistic form with definite parity and charge 
conjugation symmetry. Another important thing is how to use the relativistic wave function 
of bound state to obtain a relativistic transition amplitude, since a non-relativistic transition 
amplitude even with a relativistic wave function will lose the benefit of relativistic effects 
caused by the relativistic wave function. The Mandelstam formalism is well suited for the 
computation of relativistic transition amplitude, and we begin with this formulism to give 
a formula of the transition amplitude. 

In Sec. II, we give theoretical details for the transition amplitude in Mandelstam for- 
malism and the corresponding wave function with a well defined relativistic form. In Sec. 
Ill, the full Salpeter equation is solved, and the mass spectra and numerical value of the 
wave function are obtained. Then the two-photon decay width and full width of heavy 0~ + 
quarkonium are estimated. In Sec. Ill, short discussions and a summary are also given. 

II. THEORETICAL DETAILS 

According to the Mandelstam 10'] formalism, the relativistic transition amplitude of a 
quarkonium decaying into two photons (see figure 1) can be written as: 




4 tr I x(q) foSipi - h) h + h S(pi - k 2 ) fa 




(1) 



where k±, k 2 \ £i, e 2 are the momenta and polarization vectors of photons; e q — | for charm 
quark and e q = | for bottom quark; pi and p 2 are the momenta of constitute quark and 
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antiquark; x(q) is the quarkonium Bethe-Salpeter wave function with the total momentum 
P and relative momentum q, related by 



p t = aiP + q, aci 
p 2 = a 2 P - q, a 2 



mi 



mi + m,2 
mi + wi2 



Since p w 
mass system 



22Q. 
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M, the approximation p 10 = p 20 = ^ is a good choice for the equal 



13| . Having this approximation, we can perform the integration over q 
to reduce the expression, with the notation for the Salpeter wave function ty(q) = J -^-x{q)i 
to 



1 



n— M 



■l T tl 



(2) 



Here the relativistic Salpeter wave function ty(o) of + state with a definite parity (— ) and 
charge conjugation (+) can be written as [l^riq]: 



mi 



75 



(3) 



where u>i = \J m\ + q 2 and cj 2 = ym^ + <?■ The wave function <pi(q), ¥2(0) an d bound state 
mass M can be obtained by solving the full Salpeter equation with the constituent quark 
mass as input, and they should satisfy the normalization condition: 

dq 8u>i 



-<Pi(q)Mq) = 2M . 



(4) 



(2tt) 3 mi 

Putting wave function into the amplitude Eq. (j2J) and performing the trace, the 
amplitude becomes 
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(pi - h) 2 (Pi ~ k 2 ) 2 



M J (2tt) 3 

With this relativistic amplitude, the two photon decay width can be written as 

2 



(5) 



r(0" 



77 ) = 127ra 2 4M || j^Mo) 



+ 



1 



(pi - ki) 2 (pi - k 2 ) 



(6) 



2 

where a — |-. One can easily check that in the non- relativistic limit (by removing the 
dependence on the relative momentum q) the decay width depends on the wave function at 
the origin. 
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FIG. 1: Two-photon annihilation diagrams of the quarkonium. 

The two gluon decay width of quarkonium can be easily obtained from the two photon 
decay width, with a simple replacement in the photon decay width formula 



4 ~ 2 ~ -a, 2 . 



e q a 
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(7) 



That is: 



r((T+ - gg) = jjvr^M j J j^y^M 



1 1 

+ 



(pi - h) 2 (pi - ^) 2 



(8) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

In our previous works jlj, Q > the full Salpeter equation has been solved and the cor- 
responding eigenvalue and the wave function Eq. (jSJ) have been obtained numerically. We 
will not show the details of the calculation here, but only give the final results; interested 
readers can find them in Refs. 15^ . 

When solving the full Salpeter equation, we choose a phenomenological Cornell potential. 
There are some parameters in this potential including the constituent quark mass and one 
loop running coupling constant. In previous paper Ref. Q, the following best-fit values of 



input parameters were obtained by fitting the mass spectra for heavy meson states: 

a = e = 2.7183, a = 0.06 GeV, V = -0.60 GeV, A = 0.2 GeV 2 , A QCD = 0.26 GeV and 

m c = 1.7553 GeV. 

With this parameter set, we solve the full Salpeter equation and obtain the mass spectra 
and wave functions of quarkonium. We first fit the mass of M Vc = 2.979 GeV and then get 
the mass of r]' c : M n > = 3.566 GeV, which is a little lower than the recent experimental data 
M v > c = 3.6308 ±0.0034 ±0.0010 GeV and My = 3.654 ±0.006 ±0.008 GeV 0. With the 



obtained wave function and Eq. (jBJ), we calculate the decay width of r\ c — > 27 and rf c — > 27, 
with the result: 

rfo c -> 77) = 7.14 KeV, (9) 
r(f?c 77) = 4 - 44 K eV . (10) 

To give the numerical analysis of two-gluon decays, we need to fix the value of the 
renormalization scale \x in a s (/x). In the case of rj c we choose the charm quark mass m c as 
the energy scale and obtain the coupling constant a s (m c ) = 0.36 Q- The corresponding 
two-gluon annihilation rates of rj c and rj' c are: 

T(Vc -> <?<?) = 19.6 MeV, (11) 
r(j£-(/<7) = 12.1 MeV. (12) 

For the case of rjb and 77^, our previously derived input parameters in the potential should 
not work, because they were obtained from fitting data of heavy-light mesons, rjb and r]' b 
being heavy-heavy mesons, we change the previous scale parameters to m& = 5.13 GeV, 
Aqcd = 0.2 GeV, and other parameters are not changed. With this set of parameters, the 
mass of rjb is obtained as M Vb = 9.364 GeV, about 100 MeV lower than the T mass, and rj' b 
mass as M„i = 9.941 GeV. Now the coupling constant at the scale of bottom quark mass is 



a s (m b 



0.232. The corresponding decay widths are: 



r(r) b -> 77) = 0.384 KeV , (13) 

F{Vb 77) = °- 191 KeV > ( 14 ) 

r(r?fe -> gg) = 6.98 MeV , (15) 

r(t£ - 00) = 3.47 MeV. (16) 

In Table I we list our results with theoretical uncertainties, which are obtained by varying 
all the input parameters simultaneously within ±10% of the central values, and taking the 
largest variation of the results. In this table, we assume the total width of heavy quarkonium 
is dominated by its two-gluon decay rate, T tot ~ I^. The most recent theoretical predictions 
and experimental data are also shown in the same table. 

From the tables, we can see that our results of To!, agree well with other theoretical 



estimates of Re 
with Refs. Q, 



13. If 



.13, 



and rj^ with Refs. 
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our results of and r 2 ^, agree 



16l Il7l Il8j ; but our results of rj' c are larger than the theoretical predictions 
by others, but consistent with the recent experiment data 
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TABLE I: Recent theoretical and experimental results of two-photon decay width and total width. 



r' c KeV 
£2 



r& MeV 



KeV 



r^> KeV 
2l 



r '* KeV 
22 



Ours 
Miinz [9] 
Chao [13] 
Ebert [16] 
Fabiano [17] 
Gupta [18] 
PDG [2] 
BABAR [1] 
BES [19] 
CLEO [2ffl 

L3 [21] 
AMY [22] 
E760 [23] 



7.14±0.95 
3.50±0.40 
6-7 
5.5 
7.6±1.5 
10.94 



7.2 ± 1.2 



7.6(0.8)(2.3) 
6.9(1.7)(2.1) 
27(16)(10) 
6.71^(2.3) 



19.6±2.6 



17-23 



23.03 



16.11- 
34.3(2.3) (0.9) 
17(3.7)(7.4) 
27.0(5.8) (1.4) 



23. 9t\ 2 f 



4.44±0.48 
1.38±0.30 
2 

1.8 



12.1±1.3 



5-7 



17(8.3)(2.5) 



0.384±0.047 
0.22±0.04 
0.46 
0.35 
0.466±101 
0.46 



0.191±0.025 
0.11±0.02 
0.21 
0.15 



We comment that in this work we did not include the QCD radiative correction because 
we focus mainly on the relativistic corrections, though there is no doubt that the QCD 
correction is very important and an interesting topic. We have shown the uncertainties of 
our theoretical estimates by varying all the input parameters simultaneously within ±10% 
of the central values. It should also be pointed out that within these parameter ranges the 
uncertainty caused by the value of is very important because when we determine the 

total widths in our calculation we need the precise value of 

In summary, by solving the relativistic full Salpeter equation with a well defined 
form of the wave function, we obtain the mass spectra: M Vc = 2.979 ± 0.432 GeV, 
M v > c = 3.566 ± 0.437 GeV, M Vb = 9.364 ± 1.120 GeV and = 9.941 ± 1.112 GeV. 
With the help of Mandelstam formalism for the transition amplitude, we estimate two- 
photon decay rates: T(r] c -> 27) = 7.14 ± 0.95 KeV, T(rf c -> 27) = 4.44 ± 0.48 KeV, 
F(r) b -> 2 7 ) = 0.384 ± 0.047 KeV and F(r)' b -> 2 7 ) = 0.191 ± 0.025 KeV, and the total decay 
widths: T tat {r) c ) = 19.6 ± 2.6 MeV, T tot (r]' c ) = 12.1 ± 1.3 MeV, T tot (r] b ) = 6.98 ± 0.85 MeV 
and T tot (7] b ) = 3.47 ± 0.45 MeV. 
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